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SYMMETRY AND LINEAR STABILITY IN SERRIN’S 
OVERDETERMINED PROBLEM VIA THE STABILITY OF 
THE PARALLEL SURFACE PROBLEM 


GIULIO CIRAOLO, ROLANDO MAGNANINI, AND VINCENZO VESPRI 

Abstract. We consider the solution of the problem 

—Am = f{u) and m > 0 in m = 0 on E, 

where is a bounded domain in with boundary E of class , 

0 < T < 1, and / is a locally Lipschitz continuous non-linearity. Serrin’s 
celebrated symmetry theorem states that, if the normal derivative Ui, is 
constant on P, then D must be a ball. 

In |CMS2| . it has been conjectured that Serrin’s theorem may be 
obtained by stability in the following way: first, for a solution u prove 
the estimate 

Te -Ti < Cs [ m ] p 5 

for some constant Cs depending on 5 > 0, where Ve and ri are the radii of 
a spherical annulus containing E, E'* is a surface parallel to E at distance 
5 and sufficiently close to F, and [mJps is the Lipschitz semi-norm of u 
on E'*; secondly, if in addition Uv is constant on E, show that 

[m]ps = o(C'i) as (5—>■0’'". 

In this paper, we prove that this strategy is successful. 

As a by-product of this method, for C'^’’^-regular domains, we also 
obtain a linear stability estimate for Serrin’s symmetry result. Our result 
is optimal and greatly improves the similar logarithmic-type estimate of 
| ABR| and the Holder estimate of |CMV| that was restricted to convex 
domains. 


1. Introduction 

In this paper we establish a connection between two overdetermined bound¬ 
ary value problems, Serrin’s symmetry problem and what we call the parallel 
surface problem. As a consequence, we obtain optimal stability for the for¬ 
mer, thus significantly improving previous results ( [ABRj . [CMVj i. 

Serrin’s symmetry problem concerns solutions of elliptic partial differen¬ 
tial equations subject to both Dirichlet and Neumann boundary conditions. 
Let C be a bounded domain with sufficiently smooth boundary L 
(say 0 < r < 1). As shown in one of his seminal papers, [Se], if the 

following problem 

(1.1) Au + f{u) = 0 and m> 0 in 11, u = 0 on L, 

(1.2) Uu = a on L, 
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admits a solution for a given positive constant a, then 0 must be a ball. 
Here, / : [0, +oo) —M is a locally Lipschitz continuous function and v 
is the inward unit vector field to F. Generalizations of these result are 
innumerable and we just mention |BCN] . [BNV| . [GL| . |GNN| . [Re] , 

The parallel surface problem concerns solutions of dni with a level sur¬ 
face parallel to F, that is to say the solution u of dni) is required to be 
constant at a fixed distance from F: 

(1.3) u = b on F"^. 

Here, 

(1.4) F"^ = {x G H : fiQ(x) = J} with dn(a^) = min |a^ — 2/|, 

is a surface parallel to F, and b and 6 are positive (sufficiently small) con¬ 
stants. Under sufficient conditions on F”^, also in this case, if a solution of 
(HU and ()1.3p exists, then H must be a ball (see |MS2i IMS31 IShi IGMSll 
IGGS] i. In the sequel, we will occasionally refer to problem dni), (HU as 
the parallel surface problem. 

A condition like (II.3p was considered in |MS2[ IMS3j in connection with 
time-invariant level surfaces of a solution v of the non-linear equation 

vt = A4>{v) in H X (0, oo) 
subject to the initial and boundary conditions 

V = 1 on H X {0} and u = 0 on dVL x (0, oo); 

here, (f is a C^-smooth non-linearity with (f>{0) = 0 and (j)' bounded from 
below and above by two positive constants (hence, we are dealing with a 
non-degenerate fast-diffusion equaton). A (spatial) level surface E C H of u 
is time-invariant if v is constant on E for each fixed time t > 0. 

It is proved in |MS2j that for x € H 

rv{x,t) 

4f / - - df—> dnix)'^ as t^0~^, 

Jo 

uniformly on compact subsets of H. Hence, if E is time-invariant, then it 
has to be parallel to F at some distance 6. Also, it is not difficult to show 
that the function 

/•OO 

u(x) = / (j){v{x,t)) dt, X G H, 

_ Jo _ 

satisfies (HU, with f{u) = 1 and, being E time-invariant, u satisfies (jl.3l) . 
with F”^ = E, for some positive constant b. As a consequence, H is a ball 
and, in this situation, the time-invariant surface E turns out to be a sphere. 

Condition (|1.3I) can also be re-interpreted to give a connection to transnor¬ 
mal and isoparametric functions and surfaces. We recall that, in differential 
geometry, a function u is transnormal in H if it is a solution of the equation 

(1.5) = g{u) in H, 

for some suitably smooth function o' : M —)• (0, oo); the level surfaces of u are 
called transnormal surfaces. A transnormal function that also satisfies the 
first equation in (II.ip is called an isoparametric function and its level surfaces 
are isoparametric surfaces. Isoparametric surfaces in the euclidean space can 
only be (portions of) spheres, spherical cylinders or hyperplanes 1 |L-Cj . |Sg|); 
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a list of essential references about transnormal and isoparametric functions 
and their properties in other spaces includes |Boj . |Caj . |Mij . |Waj . 

Now, notice that the (viscosity) solution u of (|1.5I) such that u = 0 on F 
takes the form u{x) = /i(dist(x, 9^)), where h is defined by 



ds 

\/^ 


= t, 


t > 0. 


It is then clear that a solution of CH) satisfies (|1.3p if and only if n = b on 
{x E n : w{x) = h{S)}. Thus, we can claim that if a solution of (|l.ll) has 
a level surface that is also a level surface of a transnormal function, then 
must be a ball. 

Both problems (11.11) . (II.2p and ([O]) have at least one feature in 

common: the proof of symmetry relies on the method of moving planes, a 
refinement, designed by J. Serrin, of a previous idea of V.I. Aleksandrov’s 
m- The evident similarity between the two problems arouses a natural 
question: to obtain the symmetry of Cl, is condition weaker or stronger 

than (jl.3ll ? 

As noticed in |CMS2j and |CMj . condition (11.31) seems to be weaker than 
in the sense clarified hereafter. As m does not imply (11.21) . the 
latter can be seen as the limit of a sequence of conditions of type (|1.3p with 
b = bn and 6 = 6n and bn and <5,1 vanishing as n —>■ oo. As disi) does not 
imply (|1.3I) either, nonetheless the oscillation on a surface parallel to T of 
a solution of dni, doi) becomes smaller than usual, the closer the surface 
is to r. A way to quantitavely express this fact is to consider the Lipschitz 
seminorm 




sup 

x.i/er'*, 

x^y 


\u{x) - u{y)\ 
k - y\ 


that controls the oscillation of m on T*^: it is not difficult to show by a Taylor- 
expansion argument (see the proof of Theorem 15.2p that, if u E C'^’'^{Cl), 
0 < T < 1, satisfies dni), (na), then 


( 1 . 6 ) 


[u]p5 = ©(d^’*"'’') as (5 —)• 0. 


This remark suggests the possibility that Serrin’s symmetry result may 
be obtained by stability by the following strategy: (i) for the solution u of 
(HU prove that, for some constant Cs depending on 5, an estimate of type 

(1.7) Te - ri < Cs[u]y^s 


holds for any sufficiently small (5 > 0, where re and r* are the radii of a 
spherical annulus centered at some point o, {x E : r* < |x — o| < re}, 
containing T — this means that Cl is nearly a ball, if u does not oscillate too 
much on T*^; (ii) if in addition u satisfies (ffU, that is Uy is constant on T, 
show that 

[ulps = as S —)■ 0”^. 

The spherical symmetry of T then will follow by choosing 5 arbitrarily small. 

In |CMS2j . the first two authors of this paper proved an estimate of type 
(HU- Unfortunately, that estimate is not sufficient for our aims, since the 
computed constant Cs blows up exponentially as d —)• 0+. Nevertheless, in 
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the case examined in |CMj , we showed that our stategy is successful for the 
very special class of ellipses. 

In this paper we shall extend the efficacy of our strategy to the class of 
(^2,T-gmooth domains. The crucial step in this direction is Theorem 14.21 
where we considerably improve inequality dm) by showing that it holds 
with a constant Cs that is a 0(5“^) as <5 —>■ O'*". Thus, (II.6p will imply that 
rj = re, that is is a ball (see Theorem 15.2p . By a little more effort, in 
Theorem 15.31 we will prove that 

re-ri<C [uy]T, 

for some positive constant C. This inequality enhances to optimality two 
previous results, both also based on the method of moving planes. In fact, 
it improves the logarithmic stability obtained in [ABR] for C'^’'^-regular 
domains and extends the linear stability obtained for convex domains in 
[CMVj . We notice that, in our inequality the seminorm [uu\y replaces the 
deviation in the C'^(r)-norm of the function u,y from a given constant, con¬ 
sidered in |ABRj . Moreover, the inequality also improves [BNSTj [Theorem 
1.2], where a Holder-type estimate was obtained, for the case in which 
f{u) = —N, by means of integral identities. 

The outline of the proof of Theorem 14.21 will be recalled in Section 2: it 
is the same as that of |CMS2[ Theorem 4.1], that relies on ideas introduced 
in |ABRj . and the use of Harnack’s and Carleson’s (or the boundary Har- 
nack’s) inequalities. In Section 3 — the heart of this paper — by the careful 
use of refined versions of those inequalities (see [Baj . [BCNj l. we prove the 
necessary lemmas that in Section 4 allow us to obtain our optimal version 
of (II.7p . Finally, in Section 5, we present our new linear stability estimate 
for the radial symmetry in Serrin’s problem; it implies symmetry for m, 
thus, the new strategy is successful. 

2. A PATH TO stability: 

THE QUANTITATIVE METHOD OF MOVING PLANES 

In this Section we introduce some notation and we review the quantitative 
study of the method of moving planes as carried out in [ABR,j and |CMS2j 
(see also m)- 

Let H be a bounded domain in {N > 2) and T be its boundary; we 
shall denote the diameter of H by Dq. The distance dfi defined in (|1.4p is 
always Lipschitz continuous on H and of class in a neighborhood of T, 
if this is of class 0 < r < 1. In fact, under this assumption on T, for 

every x G T there are balls H C H such that x G dB; denote by r^ the 
supremum of the radii of such balls and set 

tQ = min r^;. 

We then denote by IHo the number obtained by this procedure, where instead 
the interior ball B is replaced by an exterior one. 

For <5 > 0, let be the parallel set (to F), i.e. 

(2.1) = {x G H : dn(x) > <5} . 

We know that if 0 < 5 < to, then each level surface F'^ of d^, as defined by 
dm), is of class and will be referred to as a parallel surface (to F). 
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The following notations are useful to carry out the method of moving 
planes and its quantitative version; for w E and /r E M, we set: 


TT^ = {x E : X ■ U! = ^} 
T-Lfj. = {x E : X • w > /i} 
( 2 . 2 ) Af, = Annf, 

x' = X — 2(x ■ uj — fi)uj 
{A^y = {x E : x' E A^} 


a hyperplane orthogonal to w, 
the half-space on the right of vr^, 
the right-hand cap of a set A, 
the reflection of x in vr^, 
the reflected cap in vr^. 


In the sequel, we will generally use the simplified notation A' = {A^)' every 
time in which the dependence on ^ is not important. 

Set A = sup{x • cj : X E fl}; if /i < A is close to A, the reflected cap {^y)' 
is contained in 12 (see [E])) and hence we can define the number 

(2.3) A = inf{|U : (fl/i)^ C for all /2 E (/U, A)}. 

Thus, at least one of the following two cases occurs f[5e].[FY]l: 

(51) Vl' = (12;),)' is internally tangent to cA2 at some point p' E cA2' \ -nx, 
which is the reflection in tta of a point p E dVtx \ 

(52) TTx is orthogonal to dQ at some point q E dQ n vr^,. 

In the sequel, ttx and I2 a will be referred to as the critical hyperplane and the 
critical cap (in the direction oj), respectively. Corresponding to the points p 
and q, we will also consider the points p^ = p + 6 v{p) and q^ = q + 5 v{q) 
for 0 < 5 < to; notice that T*^ = {p^ : p E T}. 

Let r E (0,ro). From now on G will denote the parallel set 


G = {x ^ U, : dci{x) > r}. 


In Section 4, we shall choose r appropriately. Also, to simplify notations, by 
P and Q we shall denote p'^ and q'^, respectively — two points on dG that 
will be frequently used. 

Now, the function w defined by 


(2.4) 

r(;(x) = u(x') — u{x), 

X E I2a) 

satisfies 


Aw + c(x) w = 0 

in I2 a, 

where for x E fix 

c(x) = 

j f{u{x’)) - f{u{x)) 

. u{x^) — u{x) 

if tt(x') 7 ^ u{x 

1 

[ 0 

if u{x') = u{x 


Notice that c(x) is bounded by the Lipschitz constant £ of / in the interval 
[0, ixiuxtt]. 

o 

By an argument introduced in [Sel Theorem 2] and improved in |BNV] 
(see also m), we can assume that re > 0 in 12;, and hence, by the strong 
maximum principle applied to the inequality Aw—c~{x) w < 0 with c~{x) = 
max[—c(x), 0], we can suppose that 


w > 0 in 12;,. 
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One ingredient in our estimates of Section 3 will be Harnack’s inequality: 
thanks to this result, for fixed a E (0,1), tc satisfies the inequality 


(2.5) 


suptc < S^a inf w, 

Bar 


for any ball Bj. C itx (see [GT[ Theorem 8.20]); the Harnack constant 
can be bounded by the power x/N + Vr£ of a constant only depending on 
N and a (see [GTj l. For instance, if c{x) = 0, by the explicit Poisson’s 
representation formula for harmonic functions, we have that 


sup w < 

Bar 


1 “h Cl 
I — a 


N 


inf w, 

Bar 


for any C (see |GTj and [DBGV] b 

Now, we review the quantitative study of the method of moving planes 
established in |GMS2j . partly based on the work in |ABRj . As already 
mentioned in the Introduction, the stability of the radial configuration for 
problem (|l.ll) . ()1.3I) is obtained in (II.7p in terms of the Lipschitz seminorm 
on parallel surfaces to F. 

For a fixed direction uj we consider the critical positions and the corre¬ 
sponding points p and q, as detailed in (SI) and (S2). As shown in |MS3| . 
the method of moving planes can be applied to G instead of Q, and the 
tangency points of cases (SI) and (S2) are P and Q, respectively. It is clear 
that if an estimate like (ILZD holds for G, then the same holds for 12, since 
the difference of the radii does not change. 

The procedure to obtain (HZD is quite delicate. For the reader’s conve¬ 
nience, we give an outline of it, in which we identify 8 salient steps. 

(i) Following the proof of [GMS2I Theorem 3.3], we show that the values 
of w{P)/ dist(P, tta), in case (SI), and of the partial derivative 
in case (S2), are bounded by some constant times [njac- 
By Harnack’s inequality, the smallness obtained in (i) at the points P 
and Q propagates to any point in G\ sufficiently far from dG\ (see 
[GMS21 Lemma 3.1]). 

By using Garleson’s inequality, the estimation obtained in the previous 
step extends to any point in the cap Gx (see |GMS2[ Lemma 3.1]), thus 
obtaining the inequality 

II^IIl°°(Ga) ^ C[u\qg- 

Here, the key remark is that C only depends on N^ r, the diameter and 
the C^-regularity of G, but does not depend on the particular direction 
oj chosen. 

The union of Gd^TLx with its reflection in vr^ defines a set A, symmet¬ 
ric in the direction w, that approximates G, since the smallness of w 
bounds that of u on dX. 

Since u is the solution of dni), then u{x) grows linearly with ^^(x), 
when X moves inside 12 from F; this implies that u can not be too small 
on dG = FL 

By using both steps (iv) and (v), we find that the distance of every 
point in dX from dG is not greater than some constant times [v\qq 
f |GMS2l Lemma 3.4]). This means that X fits well G, in the sense 


(ii) 


111 


IV 


(v) 


VI 
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that X contains the parallel set (related to G) for some positive 
(small) number a controlled by [u\gG l |CMS2l Theorem 3.5]). This 
fact is what we call a quantitative approximate symmetry of G in the 
direction u. 

(vii) An approximate center of symmetry o is then determined as the in¬ 
tersection of N mutually orthogonal critical hyperplanes. As shown in 
|ABR[ Proposition 6], the distance between o and any other critical 
hyperplane can be uniformly bounded in terms of the parameter a in 
item (vi) and hence of [u\gG- 

(viii) The point o is finally chosen as the center of the spherical annulus 
{x G : rj < |x — o| < re} and the estimate (11.71) follows from |ABR[ 
Proposition 7]. 

Based on this plan, to improve (ILZD, it is sufficient to work on the es¬ 
timates in step (ii). This will be done in Section 3, by refining our use of 
Harnack’s and Carleson’s inequalities. As a matter of fact, in |CMS2j we 
merely used a standard application of Harnack’s inequality, by constructing 
a Harnack’s chain of balls of suitably chosen fixed radius r. This strategy 
only yields an exponential dependence on of the constant in (II.7p . In 
[CMV| . we improved these estimates by chosing a chain of balls with radii 
that decay linearly when the balls approach P; however, this could be done 
only when is convex (or little more) and leads to a Holder type dependency 
on of the constant in (jl.7l) . 

In Section 3 instead, we use the following plan, that we sketch for the 
case id = ei and X = 0; p and q are the points defined in (SI) and (S2). 

We fix r = rn/4, so that G = For any 0 < <5 < r^/d, the values of 

w{x)/xi at the points p^ and P = p^ can be compared in the following way 

w{P)/Pi<G6-^wip^)/pl 

where C is a constant not depending on 5. Correspondingly, we prove that 

WxAQ) < Cg 6~^ Wx:,{q^), 

where Q = g”. Then, by exploiting steps (i) and (iii), we obtain that 

(2.6) w < CgS~^ [ri]p« on the maximal cap of Gx, 

where the constant Cg is the one obtained in step (iii) by letting r = r^/d, 
and hence it does not depends on S. 

Once this work is done, steps (iv)“(viii) can be repeated to find the im¬ 
proved approximate symmetry for the parallel set G, which clearly implies 
that for 17. We underline the fact that the dependence on 6 in (12.61) is 
optimal, as [CM| indicates. 

3. Enhanced stability estimates 

In this section, we line up the major changes needed to obtain (ILZI); they 
only concern step (ii). The following lemma will be useful in the sequel. 

Lemma 3.1. Letp and q be the points defined in (SI) and (S2), respectively. 

If B is a ball of radius r^, contained in 17 and such that p or q belong to 
dB, then the center of B must belong to 17^- 
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Proof. The assertion is trivial for case (S2). If case (SI) occurs, without 
loss of generality, we can assume that uj = ei and A = 0. Since (SI) holds, 
the reflected point p' lies on dil and cannot fall inside B, since p E dB and 

Ben. 

Thus, if c is the center of B, we have that |c —y| > to = |c —p| and hence 
|ci +pi| > |ci — pi|, which implies that ci > 0 , being pi > 0 . □ 

Our first estimate is a quantitative version of Hopf’s lemma, that will be 
useful to treat both occurrences (SI) and (S 2 ). 

Lemma 3.2. Let Br = {x E : |x| < i?} and, for p E OBr and s E 
{0,R), set p^ = p + sv{p) = (1 — s/R)p. 

Let c E L°°{Br) and suppose that w E C^{Br U {p}) 0 C‘^{Bn) satisfies 
the conditions: 

Aw + c{x)w = 0 and w > 0 in Br. 

Then, there is a constant A = A{N,R, ||c||oo)EI such that 

(3.1) u;(0) < ^s“^t(;(p®) for any sE(0,i?/2). 

Moreover, if w{p) = 0, then 

(3.2) w{0) < Awi,{p). 

Proof. We proceed as in the standard proof of Hopf’s boundary point lemma. 
Notice that w also satisfies 


Aw — c {x)w < 0 in Br, 


where c~{x) = max(—c(x), 0). Thus, the strong maximum principle implies 
that re > 0 in Br (unless tc = 0, in which case the conclusion is trivial). 

For a fixed a E (0,1) and some parameter a > 0, set 


v{x) 


7 -—— for aR < |x| < R; 

(a-“ - l)R-°‘ - I I - > 


notice that u > 0 in Br, v = Oon OBr and u = 1 on dBaR. For aR < |x| < i? 
we then compute that 


A / N — (N — 2)a — \xrc (x) , , „ 9 

Av-c-{x)v> - ) - —\x\-°^-^ 


(a-“ - 1)R- 


a‘^ — {N — 2)a — R'^\\c ' 
(a-" - l)ii- 


> 


\x 


-a-2 


Hence, we see that 


if we choose 


Av-c {x)v>0 in Br\ BaR, 


a = 


iV-2 +v^(A^-2)2+4i?2||c- 


With this choice of a, the function 


z = w — 


mm w 


1 


See the proof for its expression. 
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satisfies the inequalities; 

Az-c~{x)z<0 in Bji\BaR and z > 0 on d{Bji\BaR). 


Thus, the maximum principle gives that z >0 and hence that 


min w < 

OBaR 


w{x) 
v{x) ’ 


for X e Br\ BaR. _ 

Now, choose x = (since we want that p® G Br \ BaR, the constraint 
s/i?<l — ais needed); we thus have that 


(3.3) 


min w < 

dBaR 


o-“ - 1 

(1 - s/i?)-" - 1 


w{p’^) < 


a " — 1 w{p’^) 

K - 

a s 


where the last inequality holds for the convexity of the function t t 
Harnack’s inequality (|2.5p then yields 

— 1 wiT)^^ 

u^(0) < suptc < S^a inf W < RSja -• 

BaR ® 

Consequently, by chosing a = 1/2 and setting A = RS)i/2{‘2°‘ — I)/®, we 
readily obtain m- 

Finally, if rc(p) = 0, we readily obtain (13.21) from (13.ip and by letting s 
go to zero. □ 


The following result is crucial to treat the case (S2). 

Lemma 3.3. Set B^ = {x G Br : xi > 0} and T = {x G dB^ : xi = 0}. 
For any point q G BBrCiT and s G [0, R), define g® = q+sn{q) = (l — s/R) q. 

Let c G L°°{Br) and suppose w G C‘^{B~^) in C^{B~^ U T) satisfies the 
conditions: 

Aw + c{x)w = 0 and w > 0 in B^, w = 0 on T. 

Then, there is a constant A* = A*{N,R, ||c||oo) such that 
(3.4) u'xi(O) < A*s~^ Wxiiq’") for any s G (0,ii’/2]. 

Proof. As in Lemma 13.21 we can assume that rc > 0 in B^. 

Inequality ()3.4p will be the result of a chain of estimates: with this goal, 
we introduce the half-annulus A~^ = B^ \ B'^ and the cube 

Qa = {{xi,... ,xn) G : 0 < xi < 2cr, \xi\ < a, i = 2,... ,N}. 

For the moment, we choose 0 < p < R and 0 < u < R/y/N + 3, that is in 
such a way that C B^-, the precise value of p will be specified later. 

The first estimate of our chain is (13.5p below; in order to prove it, we 
introduce the auxiliary function 

x(x) = [|x|““ — xi for X G A+. 

It is clear that u > 0 in and u = 0 on dB'^-, also, we can choose a > 0 
so that Ax — c“(x) X > 0 in A'^ {a = {N + y/N'^ + AR?‘\\c~\\oo)/2 will do). 

We then consider the function w/v on dB'^: it is surely well-defined, 
positive and continuous in dB'j^ \ T ; also, it can be extended to be a con¬ 
tinuous function up to T n dB'^ by defining it equal to its limiting values 
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Wxx{x)/{p — R "|)forxETn dB'^. With this settings, w/v also turns 
out to be positive on the whole dB'^ since, on T n dB'^, it is positive by a 
standard application of Hopf lemma. 

These remarks tell us that the minimum of w/v on dB^ is well-defined 
and positive, and hence that the function 


z = w — rahi{w/v) v 
dB+ 


satisfies the inequalities 

Az — c~{x) z <0 in A'^ and 2 ; > 0 on dA'^. 


Thus, by the maximum principle, 2 ; > 0 on A'^ and hence 

min(t(;/i;) < w{x)/v{x) for every x E ^+. 
dB+ 


For s < R — p, we then can take x = + eei, with e > 0 so small that 

X E A'^, take the limit as e ^ 0 and, since w{q‘^) = v{q^) = 0, obtain the 
inequality 


minltc/u) < 
dB+ 




WxA(f) 

{R - s)-°‘ - R-^' 


Again, by the convexity of t i-A t ", we find that 


j^a+l 

min(t(;/i;) < - Wx^iq""), 

dB+ as 


and hence it holds that 


(3.5) 


w(x) 

mm - 

xedB+ xi 


< 


- - -^ wx, (g" 


The second estimate ()3.6p below shows that, up to a constant, the mini¬ 
mum in (j3.5p can be bounded from below by the value of w at the center of 
the cube Qa- To do this, we let y E dB'^ be a point at which the minimum 
in (13.5p is attained and set 


y = ( 0 ,y 2 , ■ ■ ■ ,yAf) and y = (cr,y 2 , • • •,^Af); 


notice that y and y coincide when yi = 0. 

The ball Ba-{y) is tangent to dB/^ n T at y and we can choose p such that 
Bfj{y) C thus, by applying Lemma [32] to Ba-{y) with p = y, = y and 
12 = ei, we obtain that 


w{y) < Aw{y)/yi if yi > 0, 


and 

w{y)<Awxiiy) if yi = 0, being u;(y) = 0. 

Thus, we have proved that 

w(x) 

w[y) < A mm -. 

xedB+ xi 

Moreover, if we also choose p such that Bp{y) C B 2 p{y) C since the 
point fjci E Bp{y), Harnack’s inequality shows that 

w{aei) < S3i/2w{y), 
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and hence we obtain that 

(3.6) w{aei) < AS^i /2 min — 

x&dB+ xi 

To conclude the proof, we use two estimates contained in |BCNj (see also 
|Ba) i■ First, after some rescaling, we can apply |BCNl Lemma 2.1] to the 
square Q 0-/2 obtain that 

(3.7) i(;(t(T/2ei) < tCi maxin, 

Qct/2 

for every t G (0,1), where C\ is the constant in [BCNl Lemma 2.1] that, 
in our case, only depends on N^ IjcUoo and R (by means of a). Thus, since 
tt)(0) = 0, taking the limit as t —>• O"'' gives that 

2C^ 

(3.8) li’xi(O) <- max in. 

^ Qa/2 

Secondly, we consider the cube and again after some rescaling, we use 
the Carleson-type estimate [BCNl Theorem 1.3] to obtain that 

(3.9) maxic < 2*? Bui(f7 ei), 

Q<t/2 

where, in our case, the constants B and q in [BCNl Eq.(1.6)] again only 
depend on N, R and IJcjjoo- Thus, by (13.8p we have that 

(3.10) inxi(O) <- w{aei). 

(7 

Therefore, by applying (13.10^ . (13.6p and (13.5p . inequality (|3.4p holds with 

A* = 2<i+^A^y2 B Cl ~ ^ , 

' a a 

where the constants p and a can be chosen as specified along the proof. □ 

For the treatment of case (SI), we must pay attention to the fact that the 
point of tangency p may be very close to ttx and the interior touching ball 
at p may not be contained in the cap. For this reason, we need the following 
lemma which gives a uniform treatment of all cases occurring when (SI) 
takes place. 

Lemma 3.4. Let ^ = ^lei with > 0 and set 

B+(0 = {x G : Jx-el <R,xi> 0}, T = {x G : \x-^\ < R, xi = O}0 

For p G dB^{^) \ T, define p^ as in Lemma, \S. SI 

Let c be essentially bounded on Bji{^) and suppose that w G C‘^{B^{fi)) n 
C^{B'I^{^) U T) satisfies 

Arc + c(x) re = 0 and w>0inB^{^), w = 0 on T. 

Then, there is a eonstant A"^ = A'^{N,R, IJcjjoo) such that 

(3.11) w{^)/^i < A^ s~^ w{p^)/pI for any sG(0,i?/2]. 

^Notice that T may be the empty set. 
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Proof. We proceed similarly to the proof of Lemma 13.31 with some modifi¬ 
cations. We shall still use the cube Qa, but we will instead consider the half 
annulus A'^ = \ 5^(0- 

Next, we change the auxiliary function v: 


v{x) = [|x-C| " -i? "]®i; 

of course n = 0 on dB'l^ ) and we can still choose a so large that v satisfies 
the inequality An — c“(x) n > 0 in A'^. Thus, the functiorU 


z = w — min(i/;/n) v 
dB+ 


is such that Az — c (x) z < 0 in A'^ and z > 0 on dA'^. By the maximum 
principle, we obtain that z > 0 on A^ and hence 


min (w/v) < w{x)/v{x) for every x G A'^. 
dB+(i) 


Again, by arguing as in the proof of Lemma 13.31 we find that 


(3.12) 


wix) — i? ") wipA 

mm - <--— 

x^dB^ (^) ^ Pi 


Now, to conclude the proof we will treat the cases > 2p and < 2p, 
separately. 

If Cl ^ ‘^Pi the ball B 2 p{f) is contained in B'^{f) and hence, by Harnack’s 
inequality, we have: 


MO < 

6 “ 2p 


< 


•^ 1/2 . 

- min w < 

2p aB+(^) 


P + P 

2p 


^1/2 


min 

x&dB+ii) 


w{x) 

Xl 


Thus, ()3.12p gives that 

^(0 / ^ ^ ^ {R/pT -1 w{p^) 

—— <{p + p) iOi/2-;;- 

Cl ' p Aa s PI 

If Cl < 2/3, we repeat the arguments of the last part of the proof of 
Lemma 13.31 with some slight modification. We take a point y € dB^ (C) 
at which the minimum in (j3.12p is attained and set y = {a,y 2 , ■ ■ ■ ,yN), 
y = (0,1/2, ■ ■ ■ We apply Lemma [3l^ to Ba-{y), with p = y, p^ = y and 

v = ei and, by inspecting the two cases yi > 0 and yi = 0, we obtain that 

wix) 

w{y) < A mm -. 

xedB+{^) Xl 


As before, we choose p such that Bp{y) C B 2 p{y) C B~^ and, since crei G 
Bp{y), by Harnack’s inequality we find that w{aei) < iOi /2 'w{y), and hence 

vji x') 

(3.13) w{aei) < S)i /2 A min -. 

xedB+{^) Xl 

Now, we apply [BCNl Theorem 1.3] to the cube Qa and obtain (|3.9I) as 
before. Moreover, again we use [BCNl Lemma 2.1] in the cube Qo-/ 2 i if P 


3As in the proof of Lemma 13.31 we observe that the function w/v can be extended 
continuously on T and w/u > 0 on 
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is sufficiently small, we have that < 2p < cr/2 and hence, applying (13.71) 
with t = 2 ^ 1 1 a gives that 

—— < Cl max in; 

Qit/2 

therefore, (|3.9p yields: 

^ < 2<iBCiw{aei). 

From (13.131) and (I3.12p . we conclude in this case, as well. The constant 
can be computed by suitably choosing p and a according to the instructions 
specified in the proof. □ 


4. Approximate symmetry 

In this section, we assist the reader to adapt the theorems obtained in 
[CMS2| in order to prove our new result of approximate symmetry for 17. 
First, we prove the analogue of |CMS2[ Theorem 3.3], that gives an estimate 
on the symmetry of fl in a fixed direction. 

Theorem 4.1. Let Q. C he a bounded domain with boundary F of class 
0 < r < 1, and set G = For a unit vector uj G let G\ and 

Ha be the maximal caps in the direetion oj for G and 11, respectively. 

Let u G be a solution of p.ll) and let w be defined by (|2.4I) . 

Then, for every 6 G (0,rn/8), we have that 

(4.1) < C(5“^[m]p5 on (a eonneeted eomponent of) Gx- 

Here, G is a eonstant depending on N, £, Do and the -regularity ofV. 

Proof. We point out that G is connected. Also, as already done before, we 
can assume that a; = ei and A = 0. 

Let p and q be the points defined in (SI) and (S2), respectively; P and Q 
are the points in dG already defined. 

In what follows, we chose to still denote by Qx and Gx the connected 
components of the maximal caps O;, and Gx that intersect B,.^^ 4 ^{P), if case 
(SI) occurs, and the connected components of Qx and Gx that intersect 
BtQ/ 4 iQ), if case (S2) occurs. 

Lemma [Q ensures that the interior ball of radius to touching dLl at p 
or q has its center in Qxl hence, P G Llx and Q G dLlx Fl ttx. We then 
apply [CMS2[ Lemma 4.2] with the following settings: Di = Gx, D2 = Llx, 
R = ro/4, and z = P, li case (SI) occurs, and z = Q, \i case (S2) occurs. 
Thus, we find that 

(4.2) w{x) < G w{P)/Pi for x G Gx, 
and 

(4.3) w{x) < G WxiiQ) for xGGx, 
respectively. Here, the constant G depends only on N,xq,£. and Dn. 
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If (SI) occurs, we apply Lemma [33] by letting R = xq,/^ and ^ = P (this 
is always possible after a translation in a direction orthogonal to ei), and 
from (14.21) we obtain that 

(4.4) w{x) <C A^5~^w{p^)/p\ for x ^ G\, 
for any <5 G (0,rQ/8). 

If (S2) occurs, we apply instead Lemma |33] (with ^ = 0 and R = r^/d) 
and (14.31) : we find that 

(4.5) w{x) < C A*6~^Wxi{q^) for x G Gx, 
for any <5 G (0,rQ/4). 

The rest of the proof runs similarly to that of [CMS21 Theorem 3.3], where 
the estimates of [CMS21 Lemma 3.2] should be replaced by (14.4p and (14.5p . 
For the reader’s convenience, we give a sketch of the proof with the usual 
settings (w = ei and A = 0). In particular, we show how to relate w{p^)/p\ 
and Wx^{q^) to [u]r 5 , which is the main argument of the proof. 

Let us assume that case (SI) occurs. If p\ > x^/2^ since p^ and its 
reflection about tta lie on T*^, then 

w{p^) = «((/)') - u{p^) < ho 
and hence we easily obtain that 

(4.6) w{p^)/p{ < 2c)or(^^ [u]r«. 

If p\ < ro/2, then \p^ — {p^)'\ < to, then every point of the segment join¬ 
ing {p^y to p^ is at a distance not greater than to from some connected 
component of T”^. The curve 7 obtained by projecting that segment on that 
component has length bounded by C \p^ — (p^)'], where C is a constant de¬ 
pending on to and the regularity of T*^ (and hence on the regularity of T 
since 6 < to/8). An application of the mean value theorem to the restriction 
of M to 7 gives that u{{p^)') — u{p^) can be estimated by the length of 7 times 
the maximum of the tangential gradient of tt on T"^. Thus, 

(4.7) w{p^) <2Cpi [u]r<5. 

Therefore, (021) and 031) yield the conclusion, if case (SI) is in force. 

Case (S2) is simpler. Since ei belongs to the tangent hyperplane to T^ at 
q^, we readily obtain (031) • Cl 

As outlined in Section [3 Theorem 14.11 completes steps (i)-(iii) and leads 
to stability bounds for the symmetry in one direction. Now, we complete 
steps (iv)-(viii). 

As described in steps (iv) and (vi), we define a symmetric open set X 
and show that G is almost equal to A. In order to do that, we need a 
priori bounds on u from below in terms of the distance function from dG, as 
specified in (v). As observed in [ABRj and |CMS2j . such a bound requires 
a positive lower bound for u^, on T, 

Uu > Co on r. 

If /(O) > 0, this is guaranteed by Hopf lemma. If /(O) < 0 instead, such 
a bound must be introduced as an assumption, as it can be realized by 
considering any (positive) multiple of the first Dirichlet eigenfunction (pi for 
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— A. In fact, for any n E N the function (pi/n satisfies (jl.ip with f{u) = 
Ai u, being Ai the first Dirichlet eigenvalue, and it is clear that, although 
{(pijn)^ —)• 0 on r as n —>■ oo, one cannot expect to derive any information 
on the shape of fl. 

The final stability result, step (viii), is obtained by defining an approxi¬ 
mate center of symmetry o as the intersection of N orthogonal hyperplanes 
as described in step (vii) (see also |CMS21 Proof of Theorem 1.1]). 

We can now conclude this section with our improved stability estimate 
on the symmetry of fl. 

Theorem 4.2. Let Q C be a bounded domain with boundary T of class 
0 < r < 1. Let u E be a solution of (II.ip . 

There exist constants 6,(7 > 0 and do E (0,rn/4) such that, if 

(4.8) [u]r«o < e, 

then there are two concentric balls and such that 

(4.9) Br^cnc Br, 
and 

(4.10) re - ri < C d~^ [u]ys , 
for any 5 E (0,(5o]. 

The constants e and C only depend on N, tn, JIq, Co, max^u and the 
-regularity o/T. 

Proof. Thanks to Theorem 14.11 we can repeat the argument of the proof of 
[CMS21 Theorem 4.2] in which we replace formula [CMS21 (3.15)] by (|4.ip . 
Hence, there exists two concentric balls Br* and Br* and two constants e 
(independent of 6) such that 

Br* C G C Br* and r* — r* < C 

if [ii]r'5 < Moreover, since 6 does not depend on S and 

lim [u]m = 0, 

we can find do E (0,ro) such that [u]p5 < e for any 5 E (0, do). To complete 
the proof, we observe that ()4.9p and (I4.10p hold with = r* -|- rn/4 and 
re = r*-\-xn/4:. □ 


5. Serrin’s problem 


In this section, we give a new proof of Serrin’s symmetry result and a 
corresponding stability estimate for spherical symmetry by using the im¬ 
proved stability inequality for the parallel surface problem (jl.ip . (11.31) . as 
just proved in Theorem 14.21 We need the following lemma. 


Lemma 5.1. Let Ll be a bounded domain with boundary T of class and 
set r = min(rn,iHn). Let u be of class in a neighborhood ofV and such 
that u = 0 onT. 

Then 


5 

1 — 5/r 



{5-t){r-t) 


[rit/z/jpi dt, 


(5.1) 




r — 5 
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for every S E [0, r). 

In particular, if 6 < r/2, we have that 

(5.2) Hr-5 < 2 J + J [iti/i/jr* • 

Proof. Let pi and p 2 be two points on F, so that pf = Pi + S v{pi), i = 1, 2, 
are points on F'^. It is clear that \p\ — P 2 I — Ifi “F 2 I — “ ^{P 2 )\ and 

hence: 

(5.3) \pi-pi\>{I-5/r)\pi-p2\. 

By applying Taylor’s formula to the values of u at pf and P 2 and taking 
the difference, we have that 

u{p\) - u{pI) = 5 [uu{pi) - Uy{p 2 )] + {S -t) [uuu{p\) - Uuuipi)] dt, 

Jo 

since u = 0 at pi and p 2 and being n{pl) = p(pi), i = 1,2. Dividing both 
sides by \p\ — P 2 \ and using (15.3p , gives that 


ujpj) 

\Pi 



d lui^(pi) - u,,(p2)l ^ 
1 - 6 /r IP 1 -P 2 I 
((5 - t)(r - t) 

r-S \p\-pi\ 


dt < 


[uu]r + 


1 — 6 /r 

Therefore, and hence (|5.2p follow at once. 


{ 6 -t){r-t) 
r — 6 


dl‘ 


□ 


We are now in position to prove both symmetry and stability for Serrin’s 
problem. Of course, stability implies symmetry, when the normal derivative 
of u is exactly constant on F. However, we prefer to present the two results 
separately. 

Theorem 5.2 (Symmetry). Let LI C he a bounded domain with bound¬ 
ary F of class 0 < T < 1. Let u E satisfy (11.11) and suppose 

that o holds with a > 0. 

Then LI is a ball. 


Proof. The assumed regularity of u implies that [uyy\Yt = 0 {P ^) as t —)• 0 
and hence, since (na) is in force, O tells us that (USD holds for some 
60 > 0. Thus, Theorem 14.21 can be applied and, by (15.2p . we have that 

fS 

Br^ C Ll C and r^ — Vi < 2C / {uyi,]Yt dt 

Jo 

for any 6 E (0, do). The behavior of as f ^ 0 then implies that the 

integral at the right-hand side can be made arbitrarily small. Therefore, 
Te = Tj, that implies that D is a ball. □ 

Theorem 5.3 (Stability). Let Ll C be a bounded domain with boundary 
F of class 0 < r < 1, and let u E C‘^^'^{Ll) he solution of (jl.ip . Let C 
be the constant in (|4.10l) . 
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There are two concentric balls Br^ and such that (BJI) holds with 
(5.4) re-ri <2C[uy]r. 

In particular, if o is in force with a > 0, then Q is a ball. 


Proof. The regularity of u and ()5.2p imply that (|4.8D holds for some Jq > 0- 
Thus, Theorem 14.21 can be applied and, by (15.2p . we have that (j4.9p holds 
with 


- r,: 


<2C < [uu]t + / [uuu]rt dt 


f 


for every 6 G (0, (5o); (15.4p then follows by letting 6 tend to 0. □ 


Remark 5.4. We notice that in Theorem 15.31 we are not assuming the small¬ 
ness of [Mj/]r to prove (lOD . 
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